arXiv: 1502.05175v2 [quant-ph] 24 Feb 2015 


Engineering adiabaticity at an avoided crossing with optimal control 

T. Chasseur,^ L. S. Theis/ Y. R. Sanders,^ D. J. Egger/ and F. K. Wilhelm^’^ 

^ Universitdt des Saarlandes, Saarbriicken, Germany 
^IQC and Dept, of Physics and Astronomy, University of Waterloo, 

200 University Ave. W, Waterloo, ON, N2L 3G1, Ganada 

We investigate ways to optimize adiabaticity and diabaticity in the Landau-Zener model with 
non-uniform sweeps. We show how diabaticity can be engineered with a pulse consisting of a linear 
sweep augmented by an oscillating term. We show that the oscillation leads to jumps in populations 
whose value can be accurately modeled using a model of multiple, photon-assisted Landau-Zener 
transitions, which generalizes work by Wubs et al. [New J. Phys. 7, 218 (2005)]. We extend the 
study on diabaticity using methods derived from optimal control. We also show how to preserve 
adiabaticity with optimal pulses at limited time, Ending a non-uniform quantum speed limit. 


I. INTRODUCTION 


The adiabatic theorem, which should be applied with 
care [1-3], states that if the time evolution of a quan¬ 
tum system is sufficiently slow, transitions between eigen¬ 
states can be neglected. It is thus a statement about 
an approximation rather than a rigorous theorem in the 
mathematical sense. Adiabatic quantum computing is 
one paradigmatic example of the usefulness of such time 
evolutions [4, 5]. Another application concerns the con¬ 
trol of quantum processor elements by frequency tuning. 
On the one hand, gate designs often rely on adiabaticity 
[6, 7]. On the other hand, due to spurious couplings — 
e.g. higher-order interactions beyond nearest neighbours 
to other parts of the chip [8] and undesired spurious res¬ 
onators [9-12] — and adiabatic following (i.e., reach a 
nonadiabatic sweep within the limited bandwidth of a 
realistic experiment), one desires a detour to diabaticity. 
Shortcuts to adiabaticity, i.e. arriving at the same final 
state as the adiabatic evolution but in a shorter time, 
have been investigated [13]. The study of adiabaticity in 
periodically driven systems has led to adiabatic Floquet 
theory [14-16]. The physics of adiabaticity is well cap¬ 
tured in the Landau-Zener (LZ) model which analytically 
describes the behavior of a system when linearly swept 
through an avoided level crossing [17, 18]. The generic 
nature of this model gives it a wide range of applica¬ 
tions, for instance in transitionless quantum driving of 
spins [19]. Considering more complex, i.e. non-linear, 
pulses has useful applications, such as Landau-Zener- 
Stiickelberg interferometry. This allows one to determine 
the parameters of the avoided crossing by quickly sweep¬ 
ing back and forth through it [20-22]. The LZ model 
also describes tunneling states in the tunneling model 
of amorphous solids. Here an atom can move between 
two adjacent positions separated by a potential barrier. 
Therefore, non adiabatic driving of tunneling states af¬ 
fects the dielectric constant in glasses [23, 24]. Evidence 
for interacting defects in glasses has been presented in 

[25] . Studies of the magnetisation of molecular magnets 

[26] involve a similar situation where a magnetic field is 
swept over many spins [27]. These dynamics also occur 
in molecular collisions [28]. 


Parallel to these developments are those in quantum 
optimal control where a control pulse is shaped to re¬ 
alize a specific time evolution [29, 30]. These methods 
were originally pioneered for nuclear magnetic resonance 
and have started to gain popularity in solid state quan¬ 
tum information devices. For instance in electron spin 
qubits to engineer gates [31], as well as in superconduct¬ 
ing qubits [32, 33] to address leakage [34, 35], frequency 
crowding [36, 37] and gate design [38, 39]. In such sys¬ 
tems the crucial entangling gates are realized with an 
anti-crossing [40]. Within the framework of qubits and 
optimal control, the LZ model has been used to study 
quantum speed limits [41-43] and how robust high fi¬ 
delity pulses are to uncertainties in the non-controllable 
part of the Hamiltonian [44]. 

For the detour to diabaticity we investigate the dy¬ 
namics of a LZ system under a linear sweep augmented 
by a fast oscillation. A similar system has been stud¬ 
ied in [45, 46]. We investigate how such a pulse can be 
used to engineer diabaticity. Optimal control methods al¬ 
low us to deepen this study as well as investigate pulses 
that keep the evolution adiabatic [13] which are crucial to 
quantum computing. The work is structured in the fol¬ 
lowing way. Section H discusses the LZ system when the 
linear sweep is augmented by an oscillating term. The 
analytics are in Sec. HA whilst the numerics are pre¬ 
sented in Sec. HB. The adiabatic pulses are discussed in 
Sec. HI. 


II. ENGINEERING DIABATICITY 

The LZ Hamiltonian is defined by ILlz = e(t)d'2/2 -|- 
A(3’a;/2 where dx and Ay are Pauli matrices in the basis 
states jo) and |1), which also are eigenstates in the ab¬ 
sence of the coupling A. A finite A mixes these states 
into new energy eigenstates referred to as instantaneous 
eigenstates and introduces an avoided level crossing lo¬ 
calized at £(t) = 0. 

Sweeping through that anti-crossing with constant ve¬ 
locity V following e(t) = vt results in the well celebrated 
transition probability between |0) at t —>• —oo and |1) 
and t ^ oo 
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PLZ = l-exp(-^) . (1) 

Therefore an linear sweep at high speed, v ^ , avoids 

leakage in our basis of |0) and |1) hence keeping the time 
evolution diabatic, whereas an infinitely slow sweep, v <C 
A^, keeps the system in the instantaneous ground (or 
excited) state at all times while that state is changing. 


A. Analytics of an Oscillation-Augmented Linear 
Sweep 


Consider a linear sweep augmented by an oscillation of 
amplitude A and frequency 

e{t) = vt + X cos + ip) . (2) 

The phase p determines the value of the oscillation when 
the linear sweep is on resonance with the anti-crossing 
and plays an important role as emphasized below. This 
form of £ can arise when investigating the dielectric con¬ 
stant of glasses by using an external electric field [23, 24]. 
Wubs et al. [45] study a very similar model where e is 
linear only and A varies sinusoidally in time. We will see 
that similar physics arises in our model and will connect 
to their work when appropriate. To study the time evo¬ 
lution produced by Eq. (2) we switch to the interaction 
picture by the transformation Uq = exp{—i(p(t)az/2} re¬ 
sulting in 

A / 0 /■* 

= y (^g-z0(t) 0 ) 

Using the Taylor expansion of sin and exp and the bi¬ 
nomial expansion of these terms are cast 

into the form 

OO 

g±i0(t) ^ 




( 3 ) 


Jm is the Bessel function of first kind. This expan¬ 
sion is similar to that done in [21, 47] and is motivated 
by coherent destruction of tunneling [48-50]. The ex¬ 
pansion in Eq. (3) allows one to identify term m with an 
jmj-photon emission/absorption process [51]. Using the 
property J-m{x) = (—l)”®Jm(a;) the Hamiltonian in the 
interaction picture is 
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For A = 0 the dynamics under a linear sweep are recov¬ 
ered. Under these conditions, the system undergoes a 
transition when t = 0 and the final population is given 
by Eq. (1). This jump is a result of the function vt'^ 
having an extrema at t = 0 and around this time the 
Hamiltonian is no longer rapidly oscillating. Thus if a 
Magnus expansion were done [52], one would see that 
the higher order commutators cannot be neglected; these 
contributions add up producing a jump. A Dyson series 
expansion also fails since as e{t) can be zero, there is no 
small parameter to expand in. 

Another way to see this is to note that Eq. (4) can be 
viewed as a sum of Landau-Zener Hamiltonians taken in 
the interaction picture where the zero bias point is shifted 
from £ = 0 to £ = mil. Indeed the terms /2+milt) 

can be recast into the form 




T*- 
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( 5 ) 


We thus expect many jumps to happen at intervals of 
mil/v which is illustrated in Fig. 1. This generalizes 
the results of Ref. [45] where only two separate jumps 
appear. 



Time [A ^] 

FIG. 1. Comparison of a linear sweep with one augmented 
by an oscillation. The time evolution was computed from 
—T /2 to T /2 with T = 100A~^ thus cutting off times without 
significant population changes. The linear speed was v = 
10A^, the oscillation frequency and amplitude were il = 50A 
and A = 120.24A respectively. The population jumps occur 
at t = —vnillv which here are integer multiples of hA”*^. 
The dotted line is the transition probability that would be 
expected of a first order Dyson series expansion. As expected, 
it does not predict the correct probability as the expansion is 
not valid since e can be smaller than A. 


1. Multi-Jump Model 

Building on the observations of the previous section we 
construct an approximate model that allows to analyti¬ 
cally compute the unitary matrix describing the jumps at 
mil/v. It consists of using the separation into multiple 
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Landau-Zener Hamiltonians and assume their applica¬ 
tion is non-overlapping in time, i.e., that the jumps are 
independent and that the sweep between the jumps only 
contributes a phase factor, which assumes A <C H. A 
time evolution that spans up to mo photon jumps is thus 
approximated by 


mo 

!/,= n 

m——mQ 



The terms in the matrices are similar to those from the 
LZ transition probability in Eq. (1) but jump m has 
Am = AJni(A/H) instead of A and an associated Stokes 
phase [21, 45] given by 


— 




TT 
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2v 


— mip 


under the assumption A'^/Av < 1. Note that negative 
Am occur, leading to an additional phase factor of tt. 
This model properly accounts for the jump heights as 
shown in Fig. 2 but is completely devoid of the oscilla¬ 
tions around these values as is to be expected given that 
each LZ jump matrix in f/j is time independent. 



Time [A-i] 


FIG. 2. Illustration of the multi-jump model. The red line 
shows the full time propagation. The horizontal lines show 
the state overlapp as predicted by the multi-jump model Uj. 
As can be seen this model accurately predicts the photon 
assisted LZ transitions. 


2. Engineering Diabaticity 


to a supression of population transfer yielding f/ = 1. 
This can be done if the time evolution for t > 0 reverses 
the time evolution of t < 0 which can be related to the 
Loschmidt echo [53]. In this section we connect this idea 
to our model and show how one has to carefully set the 
ratio A/H to a particular value. 

This can be achieved with high accuracy if the time 
during which the oscillation is present is chosen so that 
only the photon assisted jumps at m = ±1 occur. One 
has to set the phase ip to zero and choose the ratio A/H 
to the first zero of the zeroth Bessel function Jq, i.e. 
A/H ~ 2.4048. This suppresses the transition at t = 0 
and the photon assisted jumps at m = ±I cancel ea- 
chother. This follows from the Bessel function property 
J-i{x) = —Ji{x). To include the finite rise time of any 
electronics, we choose a linear ramp for the oscillation 
amplitude 

r A, if \t\ < ^ 

A(t) = I - \t\) if ^ < \t\ < ^ (7) 

L 0 otherwise 

Tg is the switching time. The height of the ramp is cho¬ 
sen to suppress the jump at to = 0, thus Jo(Ar/H) = 0 is 
needed implying Ar = 2.4048 H. The duration T — Tg over 
which the amplitude of the oscillation is held constant 
has to be long enough to include the to = ±1 transitions 
but the ramp duration plus switching time T -\-Ts has to 
be chosen so that there are no longer any oscillations in 
the pulse for jtj > 2Q,/v so as to prevent [to] > 2 jumps. 
This imposes (T — Tg)/2 > n/v and (T -|- Ts)/2 < 2^1/v. 
The time evolution of such a pulse is shown in Fig. 3. 
There is no jump at t = 0 and the to = ±1 jumps can¬ 
cel each other. As expected, this results in a pulse which 
leaves hardly any population in jO) namely less than 10“®. 
However this is phase sensitive.Indeed if ip ^ 0 then the 
Stokes phases of to = ±I are no longer different by tt and 
the reasoning presented above no longer holds. This is 
shown in Fig. 4 by the solid red line. At its worst, the 
phase can produce almost 16% error. This can be under¬ 
stood that only at phase 0, the tunneling events assisted 
by an odd number of electrons vanish automatically due 
to destructive interference of time-reversed paths. 

If even shorter pulses can be made, the phase variable 
can be rendered irrelevant. Indeed, for pulses shorter 
than H/t!, only the population jump at to = 0 could 
contribute as no photon-assisted processes are resonant. 
However, it is removed by the right choice of A/H. The 
error of a shorter pulse with T + Ts = 13.86A~^ is shown 
in Fig. 4 by the dashed blue line. 


Now the challenge of engineering diabaticity is to sweep 
through the avoided level crossing such that at the start 
and end of the evolution we stay in the same uncoupled 
state, i.e. the time evolution operator up to an irrele¬ 
vant global phase should be U{—T/2,T/2) = 1. In Ref. 
[45] it is shown how time-reversal anti-symmetry can lead 


B. Numerical optimal control 

The previous section showed that a linear sweep with 
an oscillation can reduce population jumps when going 
through the anti-crossing compared to a simple linear 
sweep. Here we show how a numerical optimization of 
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FIG. 3. Example of a pulse with v = 8A^, 0, = 50A and a 
time dependent \{t) given by Eq. (7). The ramp and switch¬ 
ing times are given by T -|- Tg = 3II/n -I- Tg = 21.55A“^. The 
ramp height satisfies Ar = 2.404811 so as to cancel the jump 
at m = 0. The final occupation of |0) is smaller than 10“®. 



Phase 9 

FIG. 4. The red solid line shows phase’s effect on the state 
overlap at T/2 of the pulse in Fig. 3. A wrongly tuned phase 
induces large errors. The dashed blue line corresponds to a 
shorter pulse where the oscillation is switched off before the 
m = ±1 jumps can occure. As expected, the error is much 
less sensitive to the phase. 


the four parameters of the pulse in Eq. (2) can reduce an 
unwanted transition when T > 2fl/v and without win¬ 
dowing. Next we show that starting from a linear sweep a 
time-sliced gradient search provides an optimal solution 
with a smaller oscillation but with a chirped frequency. 


1. Optimization of a Linear Sweep with a Single Frequency 
Oscillation Pulse 

When the short window of Eq. (7) cannot be created 
experimentally, diabaticity can still be engineered by op¬ 
timizing the four parameters (u, A, Q, ip) of Eq. (2). This 


is done with the Nelder-Mead (NM) simplex search algo¬ 
rithm [54] which is often used in optimal control [39, 55- 
57]. The metric to be minimized is ] (^(T)|0) |^. Pulses 
of arbitrarily low population transfer from jl) to jO) can 
be found with NM. Figure 5 shows a pulse with perfor¬ 
mance perfect down to 10“® where the optimization was 
stopped - in principle, machine precision can be reached. 
This pulse does not rely on suppressing the m = 0 transi¬ 
tion or time-reversal anti-symmetry. Instead, the param¬ 
eters found allow many photon assisted transitions but 
they produce a time evolution that leaves no population 
in jo) at the end of the pulse by tayloring their interfer¬ 
ence properly. Numerically we found good convergence 
regardless of the initial parameters, but typically our an¬ 
alytical solutions from the previous sections were used as 
initial guesses. Many high quality pulses having different 
parameters can be found. 



FIG. 5. Time evolution under the action of a pulse given by 
Eq. (2) where the parameters were optimized with the NM 
algorithm. The pulse duration was set to T = 200 A“^ and 
the optimization resulted in = 99.9718 A, A = 311.631 A, 
V = 9.99409A^ and p = 27r 0.387. At the end of the pulse, 
the error is of order 10~®. 


In an experiment the parameters of the pulse can differ 
from those intended. It is thus important to character¬ 
ize how robust the pulse is with respect to parameter 
fluctuations. To study how robust these pulses are we 
introduce small errors on the optimal values and study 
the decrease in fidelity. We find that these pulses are 
very sensitive to small errors in the parameters v and 
as shown in Fig. 6. The pulse is only sensitive to errors 
in the ratio between and v, not to these parameters 
separately. This can be traced back to the observation 
that the Stokes phase in Eq. (6) strongly depends on 
For parameters p and A we find that the optimal 
pulse is less sensitive and can tolerate roughly an order 
of magnitude more fluctuations than in v and Id. Snch 
errors in parameter values could be corrected by using 
the experiment to close the control loop [57]. 
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iog(i-|(v-(r)|i>p) 



- 2.0 - 1.0 0.0 1.0 2.0 


Relative Deviation in n x 10 ^ 

FIG. 6. Error in logarithmic scale as function of the error 
in the parameters (v, fl) relative to the optimum found by 
the NM algorithm for the pulse in Fig. 5. As can be seen 
their is a narrow valley where errors on v and fl do not affect 
the fiedelity too severly. However this valley is very steep 
signifying that the parameters and v need to have the 
right ratio. 


S. GRAPE pulses 

To go beyond the constraint imposed by the pulse 
shape of Eq. (2) we use the optimal control algo¬ 
rithm named Gradient Ascent Pulse Shape Engineering 
(GRAPE). It tries to optimize a fidelity 3> which is a 
functional of the pulse e{t). The pulse sequence is dis¬ 
cretized into N constant pixels Uj of time AT = N/T. 
GRAPE proceeds by iteratively updating the pixel val¬ 
ues according to the rule i-^- -I- eV^d) 

where VjXl> is the gradient of the fidelity function with 
respect to pixel j. Details of the procedure are given in 
[58]. The gradient is computed analytically [59]. For our 
situation we use the gate overlap fidelity 



-7.5 -5 -2.5 0 2.5 5 7.5 

Time [A'] 


FIG. 7. Pulse change due to GRAPE on an initial linear 
sweep of duration 9.6A~^ and v = TOrr^A^. The figure shows 
the central region of the pulse where the change is most pro¬ 
nounced. Time and pulse amplitude are respectively given in 
inverse and proportional units of the gap size A. 


III. QUANTUM SPEED LIMITS FOR AN 
ADIABATIC EVOLUTION 

This section considers the case when we wish to trans¬ 
fer the population between the two bare states by sweep¬ 
ing through the anti-crossing. Equivalently this corre¬ 
sponds to staying in the same energy branch at all times. 
This can be achieved if the time evolution is adiabatic. 
Time evolutions of this sort are important in quantum 
computing where one wishes to exchange a quanta be¬ 
tween two quantum elements such as a resonator and a 
qubit or two qubits [60]. Recently a protocol has been 
demonstrated that allows fast adiabatic two-qubit gates 
by making use of only control and optimal window 
functions [7]. For quantum computing, the population 
occupation and the phases can be important. 



( 8 ) 


where cto is the identity matrix and [/[e(t)] is the time 
evolution realized by the pulse e{t). 

When considering the LZ problem we use a linear 
sweep with speed v without coherent drive as initial 
pulse. GRAPE achieved the target error of 10“® by 
adding to the initial pulse a modulation with a time de¬ 
pendent frequency. This extra modulation is shown in 
Fig. 7. Its amplitude is much smaller than the fixed fre¬ 
quency case. The new oscillation is reminiscent of the 
Stiickelberg oscillations with decaying amplitude and in¬ 
creasing frequency. This can be viewed as a rotating 
frame version of the oscillation-augmented sweep - in¬ 
stead of sweeping the energy splitting, the frequency is 
swept accordingly. 


A. Adiabatic Pulses 


In this section we wish to realize the gate Udes. = 
through a pulse with boundary conditions e{±T/2) = 
±eo. The intial pulse for the GRAPE optimization is 


e{t) 


—A tan 


arctan 



erf (-^i) 


—erf (At) 


( 9 ) 


The parameter XT controls the width of the sweep hence 
A alone sets the speed at which the point e = 0 is 
crossed. The pulse presented in Eq. (9) is an already 
studied pulse form designed so that the sweep across 
the anti-crossing is done slowely, but further away from 
the anticrossing the sweep velocity increases [7]. Fig¬ 
ure 8 shows the phase insensitive fidelity, defined by 
‘bins = (I (0|?7]1) l^-b ] (l|f/|0) P)/2 , of the initial pulse in 
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Eq. (9) as a function of {eo,T). <f>ins only measures the 
population transfer. There are regions where the pulse 
does not perform as well as expected even when the gate 
time is long. 

Fidelity «t> 
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FIG. 8. Local phase insensitive fidelity for the initial pulse 
given by equation (9) for different gate times and start values 
for A = 0.04 GHz in all simulations. This shows that this 
start pulse fails to keep the population in the energy eigen¬ 
state whilst crossing the anticrossing, i.e. it is not adiabatic. 


[38] to prevent steep initial rise and final drop. As seen 
from the figure this empirical fit describes very well the 
behaviour observed for Tqsl- 
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FIG. 9. Fidelity of the pulses optimized with GRAPE as a 
function of duration and sweep range for A = 0.04 GHz. Once 
the gate time is sufficently large, any fidelity can be reached. 


B. Numerics 


To improve on the situation shown in Fig. 8 we use 
the GRAPE algorithm and study the dependence of the 
fidelity on gate time T. To preserve the adiabatic nature 
of the processs, the pulse is convoluted with a Gaussian to 
remove any fast oscillations. The fidelity to be optimized 
is the same as in Eq. (8) but with the target gate being (Tx 
instead of the identity. This guarantees that both phase 
and population are correct after the gate. Fig. 9 depicts ^ 

the final gate fidelity 3> in the {£o,T) parameter space £. 

where each pulse corresponds to the result of a GRAPE g 
optimization starting from Eq. (9). As can be seen all the 
regions of Fig. 8 that showed a loss of fidelity have been 
improved on. This is true as long as the pulse is long 
enough. Below a certain duration a form of quantum 
speed limit (QSL) is encountered and GRAPE can no 
longer improve the fidelity. The dependence of this QSL 
is further studied in Fig. 10. It shows the behavior of the 
QSL as function of A. The figure was created by choosing 
a gate time for which a good fidelity (> 99.99%) can be 
found for almost all values of £q. The data points (blue 
squares) are then fitted to 

rQSL(A) = io + ^^ (10) 

with c and Ag as fit parameters whereas tg is controlled 
by the fact that we have used ’’buffer pixels” in GRAPE 



FIG. 10. Study of the quantum speed limit after optimization 
of the pulse with the GRAPE algorithm. For each specific A 
a gate time is found for which a high fidelity can be achieved 
for almost all values of eg- The data points in blue are then 
fitted by the dashed line corresponding to Eq. (10). 
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IV. CONCLUSION 


We have analytically discussed the influence of a lin¬ 
ear sweep augmented by fast and strong single frequency 
oscillation on the dynamics of a two level system. We 
showed how such pulses can be used to engineer dia- 
baticity without resorting to large bandwidth control. 
Optimal control can help go beyond the analytic consid¬ 
erations producing pulses with machine precision level 
error. In the adiabatic study, improved adiabatic pulses 
were found with optimal control and their quantum speed 
limit was discussed. 
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